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We discuss the optimal setup for a low energy neutrino factory in order to achieve a 5er discovery 
of a nonzero mixing angle 613, a nonzero CP phase Sep, and the mass hierarchy. We explore parent 
muon energies in the range 5-16 GeV, and baselines in the range 500-5000 km. We present the 
results in terms of the reach in sin 2 813 , emphasizing the dependence of the optimal baseline on the 
true value of Sep ■ We show that the sensitivity of a given setup typically increases with parent muon 
energy, reaching saturation for higher energies. The saturation energy is larger for longer baselines; 
we present an estimate of this dependence. In the light of the recent indications of a large O13, we 
also determine how these preferences would change if indeed a large 6*13 is confirmed. In such a case, 
the baselines ~ 2500 km (~ 1500 km) may be expected to lead to hierarchy determination (Sep 
discovery) with the minimum exposure. 



I. INTRODUCTION 

While the indications of neutrino oscillations first came 
from the solar and atmospheric neutrino data, the obser- 
vations from terrestrial experiments have helped provide 
a firm footing to our knowledge of neutrino masses and 
mixing. The data from all the neutrino oscillation exper- 
iments have established that there are two independent 
mass squared differences lAm 2 ^ sa 2.35 x 10 -3 eV 2 and 
Awiji ~ 7.58 x 10 -5 eV 2 , as well as two large mixing 
angles sin 2 6> 23 ~ 0.42 and sin 2 12 « 0.306 (Ji]. The 
third mixing angle #13 is small: at the 3cr level we have 
an upper bound sin 2 6 13 < 0.044(0.035) Q([3). Recent 
indications of a nonzero #13 have been obtained at the 
T2K 3 and MINOS Q experiments, and now the global 
fits that incorporate these new data give a non-zero #13 
at ~ 3cr. However, the precise value of the #13 best fit 
point as well as the significance for #13 > still depend 
on assumptions on the analysis of data from reactor ex- 
periments 0. 

The immediate goals for neutrino oscillation experi- 
ments are the measurements of (i) the mixing angle #13, 
(ii) the CP violating phase Sqp, and (hi) the sign of 
Am^, also known as the mass hierarchy. These three 
quantities, along with the precision measurements of the 
already known ones, are necessary in order to complete 
our knowledge of the neutrino mass spectrum. 

Among the three quantities mentioned above, #13 is 
the most important since the determination of the other 
two depends crucially on the value of this parameter. 
If $13 = 0, the CP phase is an unphysical quantity, 
while the determination of mass hierarchy, though possi- 
ble in principle Q, becomes extremely challenging. If 
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sin 2 20i3 0.01, its measurement will be within the 
reach of accelerator experiments like T2K, MINOS or 
NO^A that use conventional hadron beams, for certain 
Sqp values. Indeed the recent results 0, 0] indicate 
that such a measurement may soon be possible. If on 
the other hand the — s- v e signals observed here are 
background or statistical fluctuations, then (9 13 may be 
even smaller. Reactor experiments, Double CHOOZ, 
Reno and Daya-Bay may probe this angle as long as 
sin 2 26»i3 > 0.033, 0.018, and 0.007 respectively @, inde- 
pendent of the value of ^cp • 

The most efficient way of determining mass hierarchy 
is the observation of the difference in Earth matter ef- 
fects for the two hierarchies, which is possible if #13 is 
sufficiently large. Future atmospheric neutrino experi- 
ments can achieve this task for sin 2 2#i3 ~ 0.04 at 95% 
C.L. 0. For T2K and NOz/A, the sensitivity to hierarchy 
is possible for sin 2 26>i 3 > 0.02 at 90% C.L., albeit only 
for limited range of values of Sep @] ■ 

The measurement of Sep is perhaps the most difficult 
of the three. Not only does it need a substantial value of 
6*i3, the value of <5cp itself needs to be sufficiently different 
from zero for a positive signal of CP violation. Planned 
supcrbeam experiments, which would use highly intense 
conventional beams, have a limited sensitivity to <5cp at 
90% C.L. and almost no sensitivity at 3a Q. 

While the technology for the conventional beams is 
well-established, the beam contamination inherent in 
such beams does not allow measurements accurate to 
more than a percent level. On the other hand, neutrino 
production from the decays of muons that are acceler- 
ated and stored in a ring ( "neutrino factory" ) , combined 
with detectors that can identify the charge of leptons 
produced from the neutrino interactions, has the poten- 
tial of measuring the quantities of interest even for much 
smaller 613 values. Even if the conventional beams suc- 
ceed in an measurement, it is important to confirm such 
a measurement with another type of source, just like the 
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measurements from solar and atmospheric experiments 
were later confirmed and established by terrestrial neu- 
trino experiments. A neutrino factory is thus a discovery 
machine in the worst-case scenario (measurements be- 
yond the reach of conventional beams) and a precision 
machine in favorable scenarios. 

In a neutrino factory, accelerated muons are allowed to 
decay in the long straight sections of a storage ring, re- 
sulting in a strong collimated beam. A fi + beam decays 
to give and v e . Oscillations of the v e to produce a 
fi~ in the detector giving the so called "wrong sign" muon 
signal, whereas the unoscillatcd produce the "right 
sign" or "same sign" muon signal of A detector with 
a charge identification capability can identify the two dif- 
ferent type of signals separately, allowing the determina- 
tion of P e)1 = P(v e -> Vp) and P w = P(D^ -> D^). A 
fj,~ beam would similarly lead to the measurements of 
P e _^ = P{v e -> Vy) and P w = P(i/ (U -> v^). While the 
sensitivity stems mainly from the wrong sign muon signal 
due to the appearance channels v e — > and D e — ¥ , the 
disappearance channels — > and — > £y, also con- 
tribute because of the large statistics available in these 
channels. The sensitivity of a neutrino factory may be 
up to sin 2 26*13 as low as ~ 10 -4 @. 

Since the appearance probability P efl depends on all 
the three hitherto unknown quantities, the channel i> e — > 
Vy in principle contains information on all of them, and 
hence it has been acclaimed as the golden channel. How- 
ever this advantage is masked by the fact that the value of 
#cp is completely unknown. This gives rise to degenerate 
solutions making the unambiguous determination of the 
oscillation parameters an uphill task. An elegant solution 
was provided by observing that at the distance ~7500 
km, the golden channel probability P e ^ becomes inde- 
pendent of the CP violating phase irrespective of energy, 
hierarchy and oscillation parameters jlOj . This baseline 
is the so called "magic" baseline. A neutrino factory is 
considered most suitable for a magic baseline experiment 
because of the large flux and the small beam background. 

Of course since the oscillation probability at the magic 
baseline is independent of the CP phase, there is no sen- 
sitivity to <5cp- Therefore though the magic baseline ex- 
periment is suitable for determination of hierarchy and 
6*i3, one has to consider other shorter baselines for the 
<5cp determination. Detailed energy-baseline optimiza- 
tion studies by the IDS-NF group propose two magne- 
tized iron neutrino detectors (MIND), one at a distance 
of 4000 km and another at the distance of 7500 km, with 
a muon energy of 25 GcV (Tlj . However this requires 
high acceleration of the muons, and one has to contend 
with the 1/r 2 fall-off of the flux. 

In the recent past many authors have investigated the 
prospect of having a neutrino factory of much lower en- 
ergy (4 -10 GeV), and hence a shorter baseline. This 
was termed as the low energy neutrino factory (LENF) 
[lH ■ The preferred detectors at these low energies are the 
magnetized totally active scintillator detectors (TASD) 
or Liquid Argon detectors that can detect muon change 



efficiently. Recently the MIND-type detectors for LENF 
were also considered [13j . Non-magnetic detectors for 
LENF has been considered in [l4[ . 

The most discussed baseline in the context of LENF 
is the DUSEL baseline of 1300 km. However recently it 
was pointed out in |l5| that the baseline ~ 2540 km has 
a special property that in the inverted hierarchy (IH), 
the golden channel probability P efl is independent of the 
CP phase around 3.3 GeV and hence it can be used for 
an efficient determination of hierarchy using the v !L — > v e 
channel in superbeams. In a subsequent Letter [161 ] . we 
showed that for the same baseline, the probability P e)Jl 
in the normal hierarchy (NH) also becomes independent 
of the CP phase at 1.9 GeV. Therefore we termed this 
as the "bimagic" baseline and the energies as the magic 
energies. We also observed that away from the magic 
energies, the probabilities still depend on Sep- Therefore 
an experiment at this baseline would be sensitive to all 
the three parameters if one uses a broadband neutrino 
beam from 1-4 GeV as can be obtained from a 5 GeV 
neutrino factory. A noteworthy point is that the distance 
2540 km which was motivated purely from physics con- 
siderations in [lU and El happens to be close to the 
Brookhaven-Homestake |17| and CERN-Pyhasalmi fl8j 
baselines 1 . 

We note that the LENF was initially motivated to 
study precision neutrino properties at large #13. There- 
fore, with the current indication that #13 is indeed large, 
an exploration of the potential of a LENF is worthwhile 
and timely. Since the neutrino-factory technology is still 
not well-established, the time scale at which these exper- 
iments will start is comparatively larger, and therefore 
its aim should be correspondingly higher - like the mea- 
surements of the above quantities of interest to 5a. This 
is particularly true for a quantity like the mass hierarchy, 
which is a binary measurement. A large #13 is also con- 
ducive for a measurement of Sep and it is likely that this 
aspect may play a decisive role in ascertaining which is 
the optimal baseline and energy. Of course if 6*13 happens 
to be smaller then it needs to be explored what is the op- 
timum baseline and energy for determination of the all 
the three unknowns - mass hierarchy, #13 and Sep- 

In general, optimization is a complex numerical prob- 
lem. The parameters involved are energy, baseline as well 
as the true values of the oscillation parameters. There is 
also the issue of the optimization with respect to the 
detector. The dependence of probabilities on the true 
values of 6*13, Sqp and Am^ is beyond the experimental 
control, so the reach of an experiment has to be assessed 
in the worst case scenario as far as the values of the mix- 
ing parameters are concerned. As we shall show further 
in this paper the lack of knowledge of true Sep makes 
it extremely difficult to zero in on a particular baseline 
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as "the optimal baseline" using the so called green-field 
approach, where one determines the optimal baseline by 
a numerical scan of the relevant parameter space. 

Optimization studies in the context of low e nerg y neu- 
trino factories have been carried out recently in [13j in the 
context of a MIND. In this paper we consider a detector 
that is TASD. Apart from this, one of the major differ- 
ences in their analysis and ours is that they present the 
sensitivity plots in terms of "fraction of <$cp'\ whereas 
our results are presented for all Sqp values in the range 
[0 — 2tt] and our plots reveal the specific range of values of 
(5cp for which a given baseline is sensitive to a particular 
quantity. 

The plan of the paper goes as follows. In Sec. [n] we dis- 
cuss the physics of oscillation probabilities giving us the 
bimagic condition and the deviations from this condition 
in the nearby baselines. In Sec. IIIII we discuss the exper- 
imental set up and the details of the numerical simula- 
tion. We present the results for the baseline optimization 
in Sec. IIVI and those for the muon energy optimization 
in Sec. |V] The dependence on the true values of Sqp 
and | Awful is shown by bands obtained by varying these 
parameters over their currently allowed ranges. The es- 
timation of an optimal muon energy, given a baseline, 
is outlined in Sec. I VII with a simple approximation. In 
Sec. IVIll we discuss the implications of a large measured 
9\3 value for the optimization. Sec. I Villi we summarize 
our results and comment on future prospects. 



II. MAGIC AND BIMAGIC BASELINES 



Certain properties of the neutrino flavor conversion 
probability P efl can be useful for an analytical under- 
standing of why certain baselines or parent muon ener- 
gies should work better than the others. While these 
have been pointed out earlier [ill, EH, HH; we expound 
on them here in detail, bringing out some of their most 
important features. 

In general the P efi oscillation probability can be writ- 
ten as 12011 



P,, 



\cos6 23 A s e ldcp +sme 23 A A \ 2 



(1) 



where As is the "solar" amplitude that depends on the 
solar parameters Am^i an d #12, and A a is the "atmo- 
spheric" amplitude which depends on Am 2 i and #13. 
From the above expression it is evident that the CP vi- 
olation in neutrino oscillation arises from the interfer- 
ence effects of these two amplitudes. In matter of con- 
stant density, the oscillation probability P Ve ^ tv can be 
expanded keeping terms upto second order in the small 



parameters a = Am^/Am^ and S13 as [2l[ 

,00 sin 2 [(1 - A)A1 9 . o„„ 9 sin 2 AA 
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where Sij = sin^y, Cy = cos 6ij. Also, 



A = 



2V2G F n e E„ 



Am 2 ! 



(2) 



(3) 



where Gf is the Fermi constant and n e is the electron 
number density. For neutrinos, the signs of A and A 
are positive for normal hierarchy and negative for in- 
verted hierarchy. A picks up an extra negative sign for 
anti- neutrinos. The last term in Eq. © corresponds to 
the interference term from which the CP dependence of 
the probability originates. This term also mixes the de- 
pendence on hierarchy and <5cp, as well as the depen- 
dence on #13 and Sqp, leading to a four- fold degeneracy 
[22^ . There is also a degeneracy between (#23,^cp) an d 
(it/2 — $23, Sqp) [Hj]. Together, this eightfold degeneracy 
makes the determination of the oscillation parameters 
ambiguous. It was noticed in [l(J that the CP depen- 
dence of the probability can be avoided if one has 



sin(AA) = 



(4) 



which corresponds to the vanishing of the solar amplitude 
As and hence the interference term. As a result P eM 
becomes independent of the CP phase Sep as well as the 
solar parameters. This condition is obeyed at the so- 
called magic baseline (L ~ 7500 km) for all E v and for 
both the hierarchies. 

However, the dependence on the CP phase also van- 
ishes when the atmospheric amplitude A a vanishes [20| . 
which corresponds to the condition 



sin[(l - A) A] = . 



(5) 



However one notes that unlike the magic baseline condi- 
tion this condition depends on energy as well as hierar- 
chy. Therefore for a particular hierarchy and a particular 
baseline one can find a set of magic energies where the CP 
dependence in the probability vanishes. If we consider 
inverted hierarchy (IH), condition for no Sep sensitivity 
(IH-NoCP) can be written as 

(1 + |i|) • |A| = nTr (6) 

with the integer n > 0. The magic energies are given as 

IH 1.27|Am 2 i|L 



E. 



magIC wr - 1.27 \K\pL ' 



(7) 



Here the quantity K has been defined such that Kp = 
2v / 2GF«e, where p = p(L) denotes the average matter 
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FIG. 1: (Color online) Graphically solving the pair of magic conditions. Left panel: E^ gic (solid) and -Emax (dashed) as a 
function of L for different n, m. Right panel: £™ ic (solid) and (dashed) as a function of L for different n, m. The 

intersection points of the curves show the baselines and the energies at which the conditions for maximum hierarchy sensitivity 
are satisfied. 



density for the baseline L. For a given baseline L, at 
these energies the probability P efl for IH is independent 
of the CP phase as well as 6*13 and only the 0(a 2 ) term 
contributes. The hierarchy dependence of the magic en- 
ergy can be utilized to maximize the hierarchy sensitivity 
by demanding that sin [(1 — A) A] = ±1 for normal hierar- 
chy (NH) at the same time. We will refer this condition 
as NH-max. With this condition, the number of events 
are enhanced due to the first term in Eq. @. Also the 
CP dependence is retained in the NH probability. The 
condition for maxima in NH can be expressed as 



(1-L4|HA| = (m-l/2)n, 
when m is any integer. This gives 



E 



XH 



l.21\&ml±\L 



(m- 1/2)tt + 1.27 \K\pL 



(8) 



(9) 



It may be conjectured that we will have maximum hier- 
archy sensitivity if 



E m ■ = 

magic max 



(10) 

From the condition in Eq. (|10[) we arrive at the baseline 

pL(km g/cc) « (n- m + 1/2) x 16300 . (11) 

Note that the relevant L in Eq. (jlip is independent of 
any oscillation parameters as in the case of the magic 
baseline. However unlike the magic baseline this condi- 
tion will be satisfied only for particular values of energy, 
given by E m ■ = E NH 

o lvVji± "-v -^magic -^max' 

Alternatively if we demand no sensitivity to CP phase 
in NH (NH-NoCP), we get 



(1-|A|).|A| = U7T 



(12) 



which gives 



E, 



NH 

magic 



l^lAmiilL 
7T + 1.27 \K\pL 



(13) 



where n is any non-zero integer. For a given baseline L, at 
these energies the probability P efi for NH is independent 
of the 5cp as well as 6*13, and P e ^ becomes 0(a 2 ). The 
condition for maxima for IH (IH-max) gives 



(l + |i|)-|A| = (m-l/2)7r 
where m is a positive integer. This gives 
l^lAm^lZ 



max 



(m- 1/2)tt- 1.27 \K\pL 



Demanding 



rlH 



E Nn . =E L 

magic max 



(14) 



(15) 



(16) 



we get the same equation for the baseline L, as in 
Eq. (fTTj) . except for an overall negative sign which can 
be attributed to the different regions of validity for n, m. 

Fig. [JJ demonstrates the existence of bimagic base- 
lines along with the corresponding magic energies, for 
n, m = 1,2,3. The left panel shows the solutions for 
IH-NoCP and NH-max, while the right panel shows NH- 
NoCP and IH-max solutions. It is clear from the fig- 
ure that the baseline L ~ 2540 km is the shortest base- 
line that satisfies both the pairs of conditions simulta- 
neously. It is not a surprise, since the equation for the 
baseline Eq. (fTTj) is similar for both the pairs of condi- 
tions. This baseline is therefore termed as bimagic flfjj ]. 
For the pair IH-NoCP and NH-max, the solution is ob- 
tained from Eq. pip with n — m and an average matter 
density of p ~ 3.2 g/cc. The corresponding magic ener- 
gies can be obtained from Eq. (O or Eq. © to be 3.3 
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FIG. 2: (Color online) Conversion probability P efl for L = 1500 km, 2500 km, 3500 km and 4000 km. The bands correspond 
to Sep € [0, 2ir]. Other parameters are as given in Eq. (fT7|) . The red (thick) line corresponds to O13 = 0. 



GeV for (n,m) = (1,1), 1.4 GeV for (n,m) = (2,2), 0.89 
GeV for (n, m) = (3,3) and so on. Clearly the magic 
energies decrease with increasing n,m. Larger values of 
n, m are less and less practical since the flux at low en- 
ergies, as well as the efficiency of detection, is typically 
lower. The solution for the pair NH-NoCP and IH-max is 
similarly obtained, the magic energies here are 1.97 GeV 
(n = 1, to = 2) and 1.09 GeV (n = 2, m = 3). 

As can be seen from the figure, the baseline of 6172 
km also satisfies the two pairs of conditions. With 
p(6172 km) = 3.955 g/cc, IH-NoCP and NH-max are 
satisfied for n — m = 1; the magic energies in the rele- 
vant range are 4.8 GeV for (n, to) = (2, 1) and 2.66 GeV 
for (n, to) = (3, 2). For NH-NoCP and IH-max, the inter- 
esting energy is 3.42 GeV, obtained with (n, m) = (1, 3). 
This baseline therefore also deserves the title "bimagic" . 
So do the longer baselines of 8950 km and 10690 km, 
which are not shown in the figure. However for the pur- 
pose of numerical optimization studies in this paper, we 
restrict ourselves to baselines in the range 500-5000 km, 
since longer baselines imply a lower flux, following the 
l/r behavior. 

Even though solving the bimagic conditions given in 
Eqs. (fTU|) and (fTTH) using the Preliminary Reference Earth 
Model (PREMTpij profile we get the exact values of the 
bimagic baselines with the magic energies, it is observed 
from Fig. [1] that the values of magic energy for one hier- 
archy and the maximum energy for the other hierarchy 
move away from each other rather slowly on cither sides 
of the intersection point, as L is varied. Moreover, cur- 



rently there is ~5% error on IAto^I at la as well as 
uncertainties associated with the density profile of earth. 
So the excellent hierarchy sensitivity of the bimagic base- 
lines, attributed to its bimagic property, is expected to 
be there even if we move slightly away from these specific 
baselines. There already exist a few possible baselines of 
similar magnitudes: (i) Brookhaven to Homcstakc dis- 
tance is exactly 2540 km, and (ii) CERN to Pyhasalmi 
(proposed site for the LENA detector) distance is 2288 
km. 

Fig-IHshows the probability P efl for different baselines, 
for sin 2 6*13 = 0.01. In this and all other plots, we have 
solved the exact neutrino propagation equation numeri- 
cally using the PREM profile. The probability P efl shown 
in Fig. [2] shows the presence of the bimagic properties for 
L ~ 2500 km clearly. Around this distance, the maxi- 
mum of NH appears at the same energy as the minimum 
of IH and vice versa, enhancing the hierarchy sensitiv- 
ity. This particular feature is absent in any of the other 
baselines shown. When one goes to a higher baseline, the 
following effects occur: 

(i) The flux decreases as l/r 2 . 

(ii) The amplitude of oscillations in NH increases while 
that in IH decreases. This would tend to increase the 
hierarchy sensitivity. 

(iii) The oscillation maxima of IH move faster to higher 
energies than that of NH, resulting in the maxima of both 
hierarchies coming closer in energies, which would tend 
to decrease the hierarchy sensitivity. 

(iv) At a given muon energy, the relevant spectral regions 
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below the highest probability peak decrease, which also 
tend to decrease hierarchy sensitivity at that energy. 
The net result is a subtle combination of all these effects. 
At the lower end of the parent muon energies E^, the 
bimagic baseline gives the best sensitivity to hierarchy, 
while at higher the optimal baseline increases. 

The figure also shows the probability for 6*13 = 0. This 
is the same for both the hierarchies since in Eq. at 
6*13=0 only the 0(a 2 ) term contributes in the leading 
order which is independent of hierarchy. The distance of 
the bands from this line gives an direct estimate of #13 
discovery potential for the corresponding baseline and 
hierarchy. The 2500 km plot in Fig. [2] suggests that #13 
discovery potential is expected to be good for NH, while 
it may not be so good for IH as the P e ^ values are lower. 
The same is true for the longer baselines. The conclusions 
get reversed for anti-neutrinos, because of the change in 
sign in A in Eq. $Z§. However for 1500 km, the difference 
between the probability values for NH and IH is not so 
large, and hence #13 discovery potential is expected to 
be similar for both the hierarchies. Of course on top 
of the probabilities, the cross section and the 1/r 2 flux 
dependence also contribute in determining the optimal 
baseline. 

The sensitivity to the CP phase is related to the widths 
of the bands that represent the variation of the CP phase. 
The widths of the P ell bands seem to increase with in- 
creasing baselines. On the other hand, the fluxes fall 
as 1/r 2 . So the optimal baseline should emerge from a 
compromise between these two opposing factors. 

In this section we have given analytic arguments to 
motivate the desirable values for the baselines and neu- 
trino energies to determine the sign of mass hierarchy, 
and detect nonzero values for #13 and <5cp- To choose the 
optimal baseline for the low energy neutrino factory, one 
will have to perform a complete numerical study, which 
we do in the next sections. 



III. DETAILS OF NUMERICAL SIMULATIONS 

As the detector, we use a 25 kt TASD with an energy 
threshold of 1 GeV. We choose a typical Neutrino factory 
setup with 5 x 10 21 useful muon decays per year, which 
is of the same order as in the setup considered in [25|, [2(| . 
We consider the running with both the polarities, each 
for 2.5 years. So we have a neutrino flux consisting of 
Dp and v e when we consider running with fi + , while it 
becomes and v e when running with negative polarity. 

We assume a muon detection efficiency of 94% for en- 
ergies above 1 GeV, 10% energy resolution for the whole 
energy range and a background level of 10 -3 for the 
v e — > and — > channels. Detection of v e or v e is 
not considered in this study, which seems to have a very 
small effect when the initial flux is as large as above [26| . 
A 2.5% normalization error and 0.01% calibration error, 
both for signal and background, have also been taken into 
account throughout this study. The detector characteris- 



tics have been simulated with GLoBES (27j . We also use 
the prescriptions for priors and marginalization in-built 
in GLoBES. 

Our main goal is to find out the optimal baseline as well 
as the optimal parent muon energy E^ for that baseline. 
In order to optimize the baseline, we choose three repre- 
sentative energies for the parent muon: 5 GeV, 7.5 GeV 
and 10 GeV, and vary the baseline in the range 500 - 
5000 km, For optimizing the parent muon energy E^ , we 
choose three representative baselines: 1500 km, 2500 km 
and 3500 km, and vary in the range 2-16 GeV. 



IV. OPTIMIZATION OF THE BASELINE 

In this section, we present the results of our baseline 
optimization for the measurement of the three "perfor- 
mance indicators" fl3l |: neutrino mass hierarchy, discov- 
ery of 6*i3 and 8cp- We have performed the analysis at 
the parent muon energies of E^ = 5, 7.5 and 10 GeV. 
The results at intermediate energies can be extrapolated 
from the results at these representative parent muon en- 
ergies. The energy optimization will be presented in the 
next section. 

The main sources of uncertainty in determining the 
reach of an experiment are the unknown values of 5cp 
and 0i3 . We therefore focus on the influence of these 
two quantities on our results, and keep the true values of 
other mixing parameters to be fixed at 

Am^ = 7.65 x 10~ 5 cV 2 , sin 2 12 = 0.3 , 
|Am^| = 2.4xlO" 3 cV 2 , sin 2 6 23 = 0.5 . (17) 

We also explore the effects of varying lAm^l in its cur- 
rent 3tr allowed range. In this section we present three 
kinds of plots: 

(i) Type-A: plots in the sin 2 #13 — L plane for fixed true 
values of lAmf.-^ and <5cp varying in the range [0, 2tt], 

(ii) Type-B: plots in the Sep — L plane for fixed true val- 
ues of sin 2 6*i3 and lAm 2 ,^, 

(iii) Type-C: plots in the sin 2 $13 — L plane for fixed true 
values of 5cp and varying (AmU in its current 3ct range. 

In the conventional plots the reach for a particular per- 
formance indicator is often given in terms of the fraction 
of Sqp values for which the determination of a quantity is 
possible. Our Type-A plots show the reach for all possible 
<$cp values, from which the reach even for the worst-case 
<$cp values may be inferred. Moreover, from the Type-B 
plots that show the reach for all Sep values, one may also 
trivially infer the fraction of 8qp for which the quantity 
may be determined. 

For reference, in the baseline optimization plots we also 
show two vertical lines, corresponding to the baselines 
of 2540 km, the "bimagic" baseline, and 1300 km, the 
baseline that is perhaps the most studied in the context 
of the LENF. 
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A. Hierarchy determination 

In order to optimize the baseline for the determina- 
tion of hierarchy, we assume NH to be the true hierarchy 
and show the reach of $13 for which the wrong hierarchy 
(IH) can be excluded to 5a, as a function of the base- 
line. Note that the determination of hierarchy is a binary 
measurement, and hence a 5a determination is absolutely 
necessary before claiming a positive identification of this 
quantity. 

Figure[3jshows the hierarchy sensitivity with the stated 
experimental setup. In the top panel we present Type-A 
plots, where true values of sin 2 #13 are plotted along the 
vertical axis. True values of all other parameters, except 
5cp, are set to values stated in Eq. (|T7l) . To generate 
the bands in this top panel, (5cp(true) is varied over the 
complete range of [0, 2ir]. For each set of chosen values of 
mixing parameters and chosen baseline, Xmin i s obtained 
by marginalizing over all parameters with the wrong hi- 
erarchy. We have taken 4% error on each of Amfi and 
#12 and 5% on #23 and jAmfil. A 2% error has also been 
considered on the earth matter profile and marginalized 
over. 

The figures in the top panel of Fig. [3] show that the 
baseline optimization depends on the actual 5cp value 
as well as the parent muon energy E^. The widths of 
the bands, which span almost an order of magnitude in 
sin 2 0i3, arc mainly due to the variation in Sep, while the 
dependence on may be discerned from the variation 
across the three plots in the panel. For each baseline, 
the upper edge of the band gives the lowest value of #13 
for which hierarchy can be determined at 5ct in the most 
conservative case, i.e. irrespective of what the true value 
of <5cp is. 

The leftmost top panel shows that for E^ = 5 GeV, the 
reach of the experiment is optimal at L ~ 2500 km for 
most of the <5cp values, however the actual reach depends 
strongly on the true 8qp- At the optimal baseline the 
reach in sin 2 #13 varies over an order of magnitude, from 
~ 10~ 4 to 10~ 3 . The maximum reach around this base- 
line is for <5cp(true) = 37r/4, which is consistent with the 
analytic estimate obtained in [l6j . For increasing parent 
muon energies, the optimal baseline typically increases. 
For longer baselines the maximum of NH and minimum 
of IH shifts to higher energies. Although the condition 
-^max = -^magic ls n0 ^ exa ctly satisfied, the broader width 
of the oscillation curve for higher baselines contributes to 
give an enhanced sensitivity. The widths of the bands in 
the top panel are rather conservative, since they include 
all true values of (5cp. 

In reality there is a unique true value of Sep, which 
controls the reach for sin 2 #13 for a given baseline. Since 
this true value is unknown in the analysis it has to be 
kept free. In order to understand what are the specific 
values of (5cp for which hierarchy can be determined at a 
particular baseline, we also present the Type-B plots in 
the middle panel of Fig. [3J plots for fixed #13 values in 
the $cp — L plane. 



The middle-panel plot for 5 GeV muon energy shows 
that if sin 2 #13 is 0.056, then it is possible to determine 
hierarchy at ha level for all values of <5cp for L > 700 km. 
For sin 2 #13 = 0.01, one needs to go beyond 1000 km if 
<5cp > 37r/2. For sin 2 $13 = 0.001, there is some sensitiv- 
ity to hierarchy for L ~ 1100 km, but only for 8qp close 
to 7r/2. Hierarchy sensitivity at 1300 km for this value 
of #13 exists only for ~ 40% of the possible <5cp values. 
The baseline that offers sensitivity to the largest range 
of <$cp values at such low #13 is 1800-2500 km. As we go 
beyond 3000 km the range of Sep values for which hier- 
archy can be determined becomes much smaller for E^ 
= 5 GeV. However if the energy is increased, then higher 
baselines can determine hierarchy for a wider range of 
<$cp values. One striking feature, seen from these plots, 
is that there is a small range of Sep values around 37r/2 
for which hierachy cannot be determined by LENF for 
sin 2 #13 = 0.001. This can be remedied by a combination 
with another experiment that is sensitive to the region 
around <5cp = 37r/2. 

Another uncertainty in the hierarchy sensitivity of a 
given experimental setup can come from IAto^I, as the 
magic energies depend on it. To illustrate this depen- 
dence, we show Type-C plots in the bottom panel of 
Fig. [3] for three fixed Jcp(true) values: i5cp(true) = 
(no CPV), ir/2 (maximum CPV), and 37r/4 (intermedi- 
ate CPV). The true value of lAm 2 ^! is varied over its 3a 
range, giving rise to bands with a finite width. Each plot 
shows the bands for the three representative E^ values, 
stated before. 

The widths of the bands in the bottom panel of Fig. [3] 
show that the dependence of hierarchy sensitivity on 
|Am|i| is much smaller compared to that on the true 
value of Sep- From the left most figure in the panel, 
we see that for no CPV, the best hierarchy sensitiv- 
ity is obtained in a rather broad baseline regime 1500 
- 3500 km, though for the worst-case values of lAm^l, 
the range around 2500 km is preferred, where the er- 
ror due to I Am 2 ;! I is also small. If CPV is maximum 
(Sep = 7r/2, central panel) or has the chosen value of 
37r/4 (right panel), the maximum hierarchy sensitivity 
is for a comparatively narrower, but still wide, baseline 
range ~ 2300 - 3500 km, near the bimagic baseline. At 
the worst-case values of lAm^l, a baseline ~ 2500 km is 
preferred for E^ = 5 GeV, while it shifts to ^3000 km 
for higher E^ values. Among the three <5cp values con- 
sidered, the sensitivity to hierarchy is observed to be the 
worst for the scenario with no CPV, however as can be 
seen from the top panel, different values of <5cp are the 
worse-case scenarios for different baselines. As the par- 
ent muon energy increases, the optimal baseline shifts 
to higher values as expected. The sensitivity is seen to 
increase when increases from 5 to 7.5 GeV, but there- 
after it seems to saturate. This feature will be discussed 
in detail in the next sections. 

Our results in this section clearly indicate that if hi- 
erarchy sensitivity is considered as the performance in- 
dicator, and if sin 2 6»i3 < 0.001 and = 5 GeV, then 
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parameters are fixed at values given in Eq. (|17[) . The plots are generated for 2.5 years of running with each muon polarity. For 
reference, we also show two vertical lines, corresponding to the baselines of 2540 km (the bimagic baseline), and 1300 km. 



the hierarchy can be determined for a larger (~ 86%) 
fraction of possible <5cp values at 2540 km, the bimagic 
baseline. This fraction increases for higher En values. 
For sin 2 #13 > 0.01, hierarchy determination is possible 
for all values of <5cp, as long as > 5 GeV and L > 1300 
km. 



B. #13 discovery 

In this section, we apply the same analysis techniques 
as the last section to the discovery potential of #13. Note 
that though a ~ 2.5a evidence for nonzero #13 has re- 
cently been claimed by experiments 0, H and the global 
fit to the neutrino mixing parameters 0, the jury is 
still out on this and it is quite possible that the value of 



#13 is much smaller. Also, since we are looking towards a 
long-term experiment, we should not be satisfied with a 
3a evidence, but should try to gauge the potential of an 
experiment for a definitive, 5a discovery. Fig.2]shows the 
5a reach for $13 discovery for three different parent muon 
energies, in terms of the three types of plots, Type-A, B 
and C, as mentioned above. 

The widths of the bands in the top panel show that 
#13 discovery potential of a baseline depends strongly 
on the true value of <5cp- Indeed for any baseline, the 
best-case and the worst-case values of Sqp make a dif- 
ference of almost an order of magnitude in the corre- 
sponding reach for sin 2 #13. Shorter baselines are seen to 
be generally preferred for #13 discovery, as compared to 
those preferred for hierarchy determination. For exam- 
ple, for E^ = 5 GeV, the optimal baseline is ~ 800—1600 
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km, while for higher parent muon energies, it moves to 
~ 1500 — 2500 km. At all these energies, it is observed 
that the sensitivity is maximum when Sep ~ 7r/2. With 
the worst-case values of <5cp, the optimal baseline stays 
near 1500 km in the whole energy range. 

In the middle panel of Fig. 2] we present the Type-B 
plots. For sin 2 9i 3 = 0.001, the discovery of #13 is possible 
for regions outside the violet (dashed) contour. Thus for 
L < 2100 km, 6*i3 can be discovered irrespective of <5cp for 
this value of sin 2 6*13 . Higher 6*13 values can be discovered 
in the whole Sep — L plane. For sin 2 #13 = 0.0001, the 
discovery of #13 is possible for the area inside the solid 
blue contours i.e. only for shorter baselines and a limited 
range of Sep- For higher baselines like L ~ 3000 km, #13 
discovery is possible only if is high and <5cp has values 
close to 7r/2 or 77r/4. In general for such small values 
of sin 2 0i3, the baseline of ^1500 km seems to have $13 
sensitivity for a wider range of Sep values. 



In the bottom panel of Fig. |H we present the Type-C 
plots for 0i3 sensitivity. It illustrates the dependence of 
the optimal baseline on the true value of <5cp : while for 
Scp = the optimal baseline seems to be ~ 1300 — 1500 
km, for S C p = tt/2 it is - 1500 - 2000 km and for 5 C p = 
37r/4, the most efficient baseline is ~ 2500 km. 



To summarize, if #13 discovery is chosen as the perfor- 
mance indicator, baselines in the region 1500 - 2000 km 
seems to work best, the actual value of the optimal base- 
line being dependent on the true value of Sep- However 
if sin 2 6»i3 > 0.001, then a LENF with E^ in the range 5 
- 10 GeV can discover $13 at 5<r for any baseline in the 
range 500-2500 km for almost all values of true i5cp • 
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C. Sep discovery 

Next we move to the Sep discovery potential in this 
section. The top panel of Fig. [5] shows the 5cr reach of 
sin 2 0i3 for <5cp discovery, given a true <5cp value. True 
values of all other parameters have been fixed to values 
mentioned in Eq. (|17p , while the true value of <5qp is cho- 
sen as 7r/4, 7r/2, 37r/4, 57r/4, 3tt/2 and 77r/4. We do not 
choose Sep (true) too close to zero or it since the discovery 
potential is expected to be low. Finally Xmin ^ s obtained 
by marginalizing over all parameters except i5cp , for each 
sin 2 8i3 (true) . Errors on different parameters are taken 
to be the same as stated in Section II V A I 

It is observed that the best-case sensitivity for Sep is 
when Sep is near tt/2 or 3ir/2. This is expected since it is 
at these values that the CP violation is maximum. Typi- 
cally for all (5cp values and all parent muon energies, the 
baselines ~ 800 — 2000 km seem to be the most efficient. 
For Sep = 37r/4, 57r/4, on the other hand, the sensitivity 



is much lower. Nevertheless, baselines of < 2000 km are 
preferred. 

The middle panel of Fig. [5] gives the sensitivity to Sep 
as a function of baseline for three fixed values of sin 2 #13 
which are 0.056, 0.01 and 0.001 and three values of energy 
Ep = 5, 7.5 and 10 GeV. These are representative values 
and the results for intermediate values of sin 2 #13 can be 
adjudged. The figure shows that for sin 2 #13 > 0.01, if 
Sqp lies in the ranges (0.3 — 0.7)7r or (1.3 — 1.7)-7r, it is 
possible to discover nonzero Sep at 5cr for all baselines 
and all energies considered here, though smaller baselines 
are a bit more efficient. For lower #13 values, only shorter 
baselines, L < 2000 km, have the possibility of detection 
of CP violation for most of the i5cp range. Longer base- 
lines need higher muon energies and specific Sep values 
near 7r/2, 37r/2 in order to achieve the task. 

The bottom panel of Fig. [5] gives the Type-C plots for 
Sep discovery for three fixed values - ir/2, 37r/4, 37r/2. 
This figure corroborates that the sensitivity to CP vi- 
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olation discovery is highest for Sep — 3tt/2. The dis- 
tance at which the best sensitivity is reached is ~1700 
km (for S C p = 7r/2) and -2100 km ( for S C p = 3vr/2). 
For S C p = 37r/4 the range 1500 - 2500 km have best 
sensitivity. In general the sensitivity is better at higher 
energies. 

To summarize, the <5cp discovery potential at various 
baselines depends on the true sin 2 #13 . For sin 2 #13 upto 
0.001, (5cp that is not too close to or 2tt can be discov- 
ered at 5ct by the LENF setup for baselines 500 - 2000 
km. In general lower baselines can access a wider range 
of Sep values. 

V. OPTIMIZING PARENT MUON ENERGY Ey 

As seen in the last section, the optimal values of the 
baselines depend on the parent muon energies in addi- 
tion to the true values of neutrino mixing parameters. In 
this section, we choose three representative values for the 
baseline - 1500, 2500 and 3500 km - and perform opti- 
mization over the parent muon energies. We use three 
different performance indicators as before: hierarchy de- 
termination, #13 discovery and <5cp discovery. The parent 
muon energy range considered here is 2 - 16 GeV, cover- 
ing the range of the proposed LENFs. We present only 
plots of Type-A, since most of the features of the energy 
dependence may be obtained through interpolation using 
the representative values Ey =5, 7.5 and 10 GeV chosen 
in the last section. 

A. Hierarchy determination 

The top panel of Fig. [5] shows the 5<r hierarchy deter- 
mination range for the three chosen benchmark baselines. 
For each chosen true value of sin 2 $13 and Ey, we obtain 
Xmin by marginalizing over all parameters assuming IH, 
while the true hierarchy is taken to be NH. All the pa- 
rameters, except Sep, are fixed to the values quoted in 
Eq. (|17|) and the band is generated by varying Sep hi the 
full range [0,27r]. Errors on different parameters have 
been taken to be the same as described in Section HV Al 

Similar to the case of baseline optimization in Fig. |3l 
the large width of the band indicates the strong depen- 
dence of hierarchy determination potential for a given ex- 
perimental setup on the true value of Sqp- As far as the 
dependence on Ey is concerned, the sensitivity increases 
with Ey at low Ey values, and saturates at a certain 
value of Ey beyond which there is virtually no change in 
the sensitivity. Let us refer to this energy as the satu- 
ration energy, E^ at . It is observed that £^ at increases 
with increasing baseline, the values being E^ at ~ 5 GeV 
for 1500 km, Ef, at ~ 7 GeV for 2500 km and Ef, at ~ 10 
GeV for 3500 km. The sensitivity reached at E^ at also 
increases with the baseline. 

The saturation behavior indicates that one will tend to 
get significantly better sensitivity with increasing Ey and 



baseline till some limit, beyond which the gain may not 
be worth the increased acceleration required. We observe 
that E^ at depends significantly on ^cp as well, however 
this dependence is not a straightforward one. The sensi- 
tivity at a baseline of 2500 or 3500 km is always better 
than that at 1500 km, sometimes by up to an order of 
magnitude, for energies beyond ~ 6 GeV. The perfor- 
mances at 2500 and 3500 km are comparable: the latter 
is marginally better, but requires higher Ey to achieve 
its full potential. 

Thus if hierarchy determination is considered as the 
performance indicator, then a LENF can determine hier- 
archy for sin 2 #13 as low as 0.002 at a baseline L ~ 2500 
km for Ey > 5GeV, for any <5cp(true). If true sin 2 # 13 
happens to be large, then shorter baselines and smaller 
energies would be sufficient. 



B. #13 discovery 

In this section we present the result of our study of $13 
discovery potential as a function of the parent muon en- 
ergy Ey, for the three benchmark baselines. The middle 
panel of Fig. [S] shows the 5a contours when <5cp (true) is 
varied in [0, 2ir] and also for seven chosen true values of 
Sep- True values of all other parameters other than the 
displayed ones have been fixed to the values mentioned 
in Eq. (fTT)) . and x 2 lin is obtained by marginalizing over 
all except $13, with the errors stated in Section HV Al 

These figures also confirm that the discovery poten- 
tial of any experimental setup depends crucially on the 
knowledge of <5cp- The sensitivity is the best when 
Sep ~ t/2, however the conservative reach for $13 may 
be up to an order of magnitude worse. Larger baselines 
require larger muon energies to get an equivalent sensi- 
tivity, however with such a larger Ey, their reach in #13 
may be better. 

The saturation behavior as in the case of hierarchy 
determination - i.e. the sensitivity increases with Ey 
till a saturation value £^ at - is observed even for 6*13 
discovery for most Sep values, prominently for baselines 
of 1500 and 2500 km, but to a smaller extent for 3500 
km. 



C. Sep discovery 

The analysis of the optimal parent muon energy for 
Sep discovery is carried out on the same lines as that in 
Sec. IIV CI The plots in the bottom panel of Fig. [5] show 
the 5er discovery potential contours as a function of Ey, 
for six chosen true values of <5cp, for each of the three 
representative baselines. 

It is observed that the best sensitivity would be reached 
when <5cp = 7r/2. While in this optimistic scenario the 
sensitivity of a 3500 km baseline with Ey > 8 GeV will be 
better, for the worst-case values of Sep a baseline of 1500 
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km is preferred, with > 5 GeV where the saturation 
is reached. 



VI. OPTIMAL PARENT MUON ENERGY FOR 
A GIVEN BASELINE 

The two design parameters of an LENF that can in 
principle be controlled are the baseline and the parent 
muon energy that determines the neutrino fluxes at the 
source. Ideally one can look for the optimal combination 
of both of these, however it may not always be practical. 
In particular, considerations behind choosing a baseline 
involve factors like the location suitable for an acceler- 
ator and a location where a undergraound laboratory 
for neutrino detection can be built. Apart from scien- 
tific considerations, this involves geography, economics 
as well as sociology. A survey of such pairs of locations 
was recently carried out [l3| . Given a baseline, the choice 



of the parent muon energy is relatively straightforward, 
limited mainly by technological considerations. The en- 
ergy of the muon beam can then be chosen based on the 
optimality analysis. 

As has been noticed multiple times in the previous two 
sections, the performance at a baseline typically increases 
with increasing parent muon energy (in the LENF 
range), till it saturates at a particular value £^ at . Given 
that increasing the muon energy is associated with addi- 
tional costs, a desirable thing to do would be to choose 
the value of E^ at or near the value of E^ at . In this 
section, we obtain the E^ at values for arbitrary baselines 
using a simple approximation. It turns out that this ap- 
proximation matches the numerical £^j at values observed 
in the previous section. 

We expect that the main factor influencing the effi- 
ciency at a given energy will be the number of wrong-sign 
muon events. While the actual numbers will depend on 
the detector characteristics, we estimate this efficiency 



13 



> 
<u 

o 

— 




500 



1500 



2500 3500 

L(km) 



4500 



FIG. 7: (Color online) The optimal value E° p as a function 
of the baseline, for normal hierarchy and with sin 2 813 = 0.001 
and 0.01. 



simply through the "quality factor" 

/ ^^P efl a„^dE^ 



Q 



f$ Ve a^dE Ve 



(18) 



where 3>i, e is the flux of v e at the source and a v is 
the total charged-current cross section of at the de- 
tector. We have normalized the events to a complete 
v e — > conversion. Since the parent muon energy at 
the source determines the neutrino spectra at the source 
completely, and the baseline determines the oscillation 
probability completely (modulo our knowledge about the 
mixing parameters), the quality factor is known for a 
energy-baseline combination. The optimal energy E° pt is 
the one for which this quality factor is maximum. There 
will be a spread in E° pt due to the uncertainties in mixing 
parameters, mainly #13, lAm^l and Sep- 

Fig. [7] shows the E° pt values as function of baseline L, 
for sin 2 0i3 = 0.001, 0.01. The width is obtained due to 
variation of E° pt as [Amf-J is varied in the current 3<r 
range and Sqp in [0, 2ir}. This gives a good estimation of 
what E^ values will be optimal corresponding to a given 
baseline, and can work as a rough guideline for determin- 
ing the optimal range of E^ , once the baseline L has been 
determined from other considerations. Then the upper 
limit can be used to estimate £^ at for a given baseline. 
Going beyond -E^ at would not improve the performance 
of the setup. 



VII. IF 0i3 IS LARGE 

The indications of a large #13 value from the exper- 
iments [H, |i[ and the global fit to neutrino data [2, [|[ 
is a good news for the measurements of other neutrino 
parameters too. In particular, it makes it easier to de- 
termine the mass hierarchy, and makes the measurement 
of Sep possible. 

The figures presented in the previous section already 
give us some idea about what can be the optimal ener- 
gies and baseline for hierarchy and Sep for sin 2 #13 in the 



above range. 



If sin 2 6*i3 is close to its present best-fit, then the LENF 
set up considered in this paper can determine hierarchy 
irrespective of what is the true 5cp for baselines <; 1000 
km, with a 2.5 years exposure with each of the polarities 
and 5 x 10 21 useful muon decays, as shown in Figure [3] 
If true Sep is not too close to 37r/2 then even baselines ^ 
700 km would be sufficient for this purpose. The task of 
Sqp determination will also be easier: Sqp in the range 
(0.2 — 0.8)7r and (1.2 — 1.8)7T can be determined at ha for 
baselines of 500-2000 km for E^ in the range 5-7 GeV, 
as can be seen from Figure [5] 

We can then be more ambitious and try to optimize 
the setup in order to get a ha determination of hierarchy 
and Sep with the minimum exposure. This is shown in 
Fig. [SI where we take the range of sin 2 6*13 to be [0.013, 
0.028], the ler range given in [l|, in the prior and then 
marginalized over #13. 

It is clear from the top panel of the Fig. [H] that for 
any parent muon energy, the baseline L ~ 2500 km is 
the optimal one for the determination of mass hierarchy. 
As mentioned before, this is near the bimagic baseline 
indicated in [l(|. For a baseline of 1500 or 3500 km, the 
exposure needed may be an order of magnitude larger, 
depending on the Sqp value. 

The bottom panel of Fig. [8] shows that for any E^ 
a baseline ^1300 km will show the best sensitivity for 
Sqp discovery, for most of the true values of Sep- How- 
ever, the "bimagic baseline" 2540 km will have compa- 
rable sensitivity with an exposure ~1.5 times larger. If 
i5cp(true) happens to be close to 57r/4, ^2500 km may 
show a better sesitivity compared to 1300 km. Note that 
the exposure needed for i5cp discovery at ha is typically 
one order of magnitude higher than that for hierarchy 
determination. 



VIII. SUMMARY AND CONCLUDING 
REMARKS 

The neutrino mixing angle $13, the CP-violating phase 
Sep and the hierarchy of neutrino mass eigenstates are 
the three quantities whose measurements are crucial in 
order to complete and confirm our current picture of neu- 
trino mixing and oscillations. The next generation of 
long-baseline neutrino experiments therefore justifiably 
consider these measurements as their primary goals. 

In this article, we go beyond the conventional beam ex- 
periments and try to optimize the setup for a low-energy 
neutrino factory (LENF) - where the energy of the parent 
muon is less than ~ 15 GeV, and consider a magnetized 
totally-active scintillator detector (TASD) that can iden- 
tify muon charge, for these measurements. In the worst 
case scenario i.e #13 turning out to be smaller than the 
reach of the current reactor and superbeam experiments 
such setups can be used as discovery machines. On the 
other hand if #13 is confirmed to be in the currently indi- 
cated range then one can use such machines for precision 
studies with the aim to measure the three quantitites at 



14 



a 

1 1 
o 

3 



NH 

E^ = S GeV 



Hit' rare hv 
determination 



<™ (T2K * MINOS result) 

2 o.i 



71/4 
71/2 

3t^4 
5JI/4 



nH 

K M = 5 GeV 




NH 

E M = 7.5 GeV 



, m*4S 

f Hierarchy ] S --^ - 

[determination J ( T2K + MINOS result) 



\ \ 

* \ \ \ 


NH 




E = 10 CcV 




^Detccto 


= 25kl 


( Hierarchy **S^ 


W 


1 determination J 


(T2K + MINOS result) 



500 1500 2500 3500 4500 500 
L (km) 
.100- rr — . ' ..- 

>j •.. [Oqp discovery I 



1500 2500 3500 4500 500 
L (km) 



1500 2500 3500 4500 
L (km) 



V\( 8 CP 


discovery] 


NH 






7.5 GeV 






^Detector = 25 kt 






(T2K + MINOS result) 





discovery] 


NH 






E M = Ifl GeV 






M Detector = 25 kt 






(T2K + MINOS result) 



1500 2500 3500 4500 500 
L (km) 



1500 2500 3500 4500 500 
L (km) 



1500 2500 3500 4500 
L (km) 



FIG. 8: (Color online) Exposure required for 5a hierarchy determination (top panel) and Sep discovery (bottom panel), if #13 
is large. We have taken sin 2 #13 £ [0.013, 0.028] [l[. NH is taken as the true hierarchy, and all parameters except <5qp has been 
fixed at the values stated in Eq. <|17|) . 



the 5<7 level. In this paper we have focussed on the dis- 
covery potential. 

We note that some of the parameters that determine 
the efficiency of a neutrino factory are beyond our con- 
trol: the probability P e/J that controls this efficiency de- 
pends crucially on #13 and Sep- However there are two 
parameters in our control: the baseline and the energy of 
the parent muon. It is with respect to these parameters 
that we perform our optimization. 

It is of course obvious that the value of 813 would de- 
termine how efficient an experiment is. And that is true 
more or less in a straightforward way. In the absence of 
matter effects, P efJ _ is proportional to # 2 3 . In the pres- 
ence of matter effects, #13 affects the probability linearly, 
thus influencing the identification of mass hierarchy. The 
term in the probability that involves Sep is also linear in 
813. Thus at shorter baselines (low matter effects), the 
value of 0i3 affects the measurement of 6*13 itself through 
a quadratic term, while at larger baselines, it affects the 
probability through a term that contains the product of 
itself with the matter effects. 

In addition to #13, we find that even the true value of 
Sep can have a large impact on the reach of an LENF 
experiment. Indeed, the reach can change by up to an 
order of magnitude or more depending on the true value 
of 5cp- The actual value of |Am§i| also affects the effi- 
ciency, albeit to a smaller extent. Given that the actual 
value of Sep is crucial, we present our results in a way 
that bring out the impact of this parameter. We present 
three types of plots: (i) Type-A: plots in the sin 2 6*13 — L 
plane for fixed values of jAm^j and Sep varying in the 



range [0, 2tt], (ii) Type-B: plots in the Sep — L plane for 
fixed values of sin 2 813 and IAto^I, (iii) Type-C: plots in 
the sin 2 #13 — L plane for fixed values of Sep and varying 
IAto^I in its current 3ct allowed range. 

Most of the earlier work with the optimization of LENF 
presented the results in terms of "fraction of <5cp" where 
the experiment was successful. However knowing the ex- 
act range of i5cp where a certain measurement is possible 
is crucial, especially when one is thinking about combin- 
ing results from two complementary experiments: these 
experiments should be sensitive to complementary ranges 
of 5 C p- 

Our detailed observations may be found in the main 
body of this text. We would like to attract the reader's 
attention to some salient features. The following obser- 
vations refer to 2.5 years of running with each muon po- 
larity, with 5 x 10 21 useful muon decays per year. 

(i) For sin 2 6*13 > 10~ 2 , mass hierarchy can be determined 
to 5cr at almost all the baselines > 1000 km. For lower 
#13, baselines ~ 2500 km achieve the task for the largest 
fraction of (5cp, however even they may fall short when 
Sep is near 37r/2. 

(ii) Any baseline < 2000 km can discover 6*13 if sin 2 6*13 > 
10~ 3 . For smaller #13, the choice of optimal baseline de- 
pend on the actual value of Sep- However for majority 
of values of Sep the baseline range 1300 - 2500 km can 
be termed as optimal. 

(iii) The discovery of Sep is naturally harder when its 
value is near or ir. However if sin 2 6»i 3 > 0.01, the dis- 
covery is possible for a wide range of Sep values centered 
at 7r/2 and 37r/2. The range typically decreases with the 
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increase of baseline. 

(iv) At a given baseline, when the parent muon energy is 
increased, the performance typically increases up to some 
energy and then remains the same. We term this as the 
saturation behavior. 

As can be gathered from above, there does not ex- 
ist a unique "optimal" baseline or muon energy for all 
performance indicators. The baseline determination de- 
pends, in addition to scientific merit, also on geography, 
economics and sociology. Once that is determined, the 
optimization of muon energy involves mainly scientific 
and technical considerations. The saturation behavior 
mentioned above indicates that, given a baseline, it is 
preferable to have the parent muon energy as close to 
the saturation energy as possible, in order to avoid in- 
creasing the parent muon energy unnecessarily. While 
the saturation energy may be determined through a de- 
tailed simulation for a given baseline (as we have done for 
three benchmark baselines in this article), we have tried 
to come up with a simple criterion that can motivate this 
behavior and determine the saturation energy. We con- 
jecture that the most efficient energy for a given baseline 
simply depends on the number of wrong-sign muons at 
the detector. This gives us a range of optimal muon ener- 
gies for a given baseline, the range depending on the ac- 
tual values of mixing parameters. This conjucture is vin- 
dicated a posteriori by the saturation energies obtained 
at the benchmark baselines. This simple-minded anal- 
ysis gives an intuitive understanding of the relationship 
between baseline and the optimal energy. 

With the recent indications of a large value of 6*i 3 , the 
5cr measurement of this quantity may already be within 
our grasp with the current experiments. Then the neu- 
trino factory experiments have to aim higher, and the 
question to ask is what kind of setup will give us the 5a 
determination of mass hierarchy and CP violation with 
the least amount of exposure needed. Our analysis indi- 
cates that for the mass hierarchy, the baseline of ~ 2500 



km would perform the best at any parent muon energy. 
On the other hand, the measurement of CP violation is 
the most efficient around a baseline of ~ 1500 km at all 
parent muon energies, for most of the <5cp values. 

A few comments are in order. We have presented all 
our results assuming the actual hierarchy to be NH. How- 
ever note that the probabilities obtained with NH and //+ 
beam are identical with the probabilities with IH and fx" 
beam, and we have taken beams of both polarities with 
equal exposure. In LENF, both the polarities of muons 
can be accelerated at the same time, giving alternate 
bunches of /i + and fi~ . The cross sections of and 
on the nucleons are also virtually identical at the relevant 
energies, and hence our results are valid even for IH. 

Another comment is about the measurement of <5cp it- 
self. In this article, we have only considered the discovery 
potential for CP violation, i.e. we are interested in find- 
ing a nonzero <5cp- This task is, naturally, hard for small 
<5cp and impossible when <5cp=0 or n. However with 
large #13, a measurement of Sep may be possible even if 
its value is close to the CP-conserving limit. If indeed 
the value of #13 is confirmed to be large, one can aim to 
answer more detailed questions like the value of Sqp, the 
octant of 823, or the deviation of #23 from maximality. 
These quantities will not be analyzed in this work. 

Finally, the most important quantity is the one we have 
no control over: the value of #13. As far as neutrino fac- 
tories are concerned, this quantity will determine if they 
are going to be discovery machines or precision machines. 
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